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We present a microscopic calculation of the nondissipative current through a superconducting 
quantum point contact coupled to a mechanical oscillator. Using the nonequihbrium Keldysh Green 
function approach, we determine the current-phase relation. The latter shows that at certain phases, 
the current is sharply suppressed. These dips in the current-phase relation provide information about 
the oscillating frequency and coupling strength of the mechanical oscillator. We also present an effec- 
tive two-level model from which we obtain analytical expressions describing the position and width 
of the dips. Our findings are of relevance for nanomechanical resonators based on superconducting 
materials. 
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I. INTRODUCTION 

The observation of a small, albeit macroscopic, me- 
chanical oscillator in the ground state [l[ constitutes 
a milestone on the road leading to the experimental 
verification of quantum decoherence in mechanical sys- 
tems. Detection and manipulation of such nanomechan- 
ical oscillators are still a challenging issue and several 
approaches are being pursued by experimentalists. A 
nonexhaustive list of devices includes electromagnetic 
cavities in the microwave range [3, l3[ , superconducting 
qubits H, optical cavities (a-Q]' single-electron transis- 
tors |8|, |9|, and tunnel junctions [lO[. But in order to 
tackle fundamental questions of decoherence in macro- 
scopic systems, |ll| still further improvements are neces- 
sary, thus motivating the investigation of new directions 
for the detection of nanoelectromechanical systems. 

In Ref.[lO| the modulation of the tunneling quantum 
amplitude in an atomic point contact has been exploited 
to detect the mechanical fluctuations of a doubly clamped 
beam. The current through the point contact is modu- 
lated by the change in the distance between the oscillat- 
ing beam and a fixed reference metal. It has been shown 
that this kind of detector can reach the quantum limit 
of displacement detection [12, llJ] and can be allowed, in 
the experiment of Ref.[l^, to measure with a good ac- 
curacy the resonating frequency of the oscillator and its 
Brownian motion. 

It seems feasible to reproduce a similar experiment 
with superconducting leads instead of normal metal ones. 
In this case the atomic point contact forms a Josephson 
junction between the mobile and fixed leads. Joseph- 
son current in similar tunnel junctions has been demon- 
strated by using a superconducting scanning tunneling 
microscope (STM) tip [Tj]. Also the current phase rela- 
tion has been measured in Josephson junctions consisting 



of atomic point contacts [la, \T^ and carbon nanotubes 
[iTj. An alternative system with a similar modulation 
of the tunneling amplitude induced by a mechanical dis- 
placement is a suspended carbon nanotube contacted be- 
tween two superconductors in the Fabry-Perot regime. 

From the theoretical point of view, the effect of the 
modulation of the tunneling matrix clement due to me- 
chanical oscillations in Josephson junctions has been con- 
sidered in the literature, but only in the adiabatic limit of 
slow variation of the tunneling amplitude [1^ [l^ • This 
limit is adapted to the tunneling case where the only rele- 
vant energy scale is the superconducting gap A. Thus for 
oscillating frequencies much smaller than A/fi, (typically 
of the order of tens of GHz) the time dependence of the 
oscillator can be treated adiabatically with a correction 
pro portional to the time derivative of the displacement 

^- . . . . 

The situation is different for an atomic point contact 

consisting of few conducting channels, some of them with 
high transmission. It is well know [20, [21[ that in such 
junctions the supercurrent is controlled by the occupa- 
tion of the Andreev bound states (ABSs) with ener- 
gies depending on the transmission and phase difference 
across the contact. In the case of a unique conducting 
channel, there are only two Andreev states with energy 



± wa(0) = ±AJl - Tsin^(0/2) , 



(1) 



where </> is the phase difference between the two super- 
conductors, and T (0 < T < 1) the transmission coeffi- 
cient of the channel. For r ^ 1 (tunneling limit) the 
spectrum of the two-level system as a function of the 
phase difference is almost flat with a level spacing be- 
tween the channels equal to 2A. In the opposite limit, 
for T = 1, the level spacing between the ABSs depends on 
(j) and shows a minimum at </> = tt where the energy split- 
ting vanishes. Thus, the energy level splitting 2uja spans 



the range 2A > 2u!a > 2Av'l — t and may be equal to 
the mechanical resonating frequency. Moreover, for large 
enough phases the two Andreev levels can be deep inside 
the superconducting gap and an effective two-level model 
description for the contact that neglects the continuum 
part of the spectrum can be used |22l | . 

A two-level system is the simplest example of a quan- 
tum detector [23| . By tuning the energy splitting between 
the two levels to a radial frequency ujq one can measure 
the transition rates between the two levels induced by the 
coupling to an external system and thereby determine 
the fluctuation spectrum at ojq. For the two- level sys- 
tem formed by the Andreev states it has been predicted 
very recently [2J, |25| that under microwave irradiation 
of radial frequency ojq the current-phase relation of the 
Josephson junction show dips at values of cj) that are so- 
lutions of the equation Hloq = 2ujj^{(j)). At these values 
of the phase the electromagnetic field induces transitions 
between the two Andreev levels and the resulting current 
vanishes at resonance. The question thus naturally arises 
if a similar mechanism can be used to detect mechanical 
oscillations. 

In this paper we explore the effect of fast modulation 
of the transparency of a single channel quantum point 
contact on the current-phase relation (CPR). The origin 
of that modulation is the vibration of the contact, with 
oscillations in the range of hundreds of MHz. We show 
that if the frequency of the mechanical oscillator is of the 
order of the spacing between the Andreev levels, then 
the CPR differs drastically from the one obtained in the 
adiabatic case. In analogy with the microwave irradia- 
tion, we find that the modulation of the tunnel amplitude 
leads also to the appearance of dips in the current-phase 
relation that, for typical experimental situations, could 
be extremely sharp. Measurement of the current phase 
relation could thus be used to detect the resonating fre- 
quency of the mechanical oscillator. In order to describe 
this effect it is clear that one has to go beyond the adia- 
batic approach of Ref. [l^ in order to take into account 
the transitions between ABSs. 

The plan of the paper is the following. In Sec. [Ill we 
first introduce a microscopic Hamiltonian to describe the 
superconducting point contact in the presence of elec- 
tromechanical vibrations. From this model Hamiltonian, 
we compute the supercurrent through the junction in two 
ways. In a first approach. Sec. IIII Al we use an effective 
two-level model similar to the one derived in Ref. j22| . 
Within this model, we derive analytical expressions for 
the current and for the occupation of the Andreev states. 
Second, in Sec. IIIIBl we calculate the supercurrent with 
the help of the Keldysh Green functions and check the 
validity of the two-level model. In particular we present 
numerical solutions for the current phase relation. In Sec. 
lIVI wc discuss the use of this method to detect mechanical 
oscillations and present the conclusions in Sec. |Vl 



II. MODEL 

We consider a model for a Josephson junction where 
the normal state electronic hopping term is modulated 
in time at a given frequency. An example of the sys- 
tem considered is depicted in Fig. [1] It consists of two 



v{l + acos{wot)} 
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FIG. 1: Schematic representation of a Josephson junction. 
We assume that the distance between the two leads is mod- 
ulated periodically at a frequency cjo- This induces a weak 
modulation of the quantum hopping matrix element v{t). 

bulk superconductors, characterized by a BCS gap A, 
connected by a junction whose dimensions are assumed 
much smaller than the superconducting coherence length. 
This model can describe, for example, the contact be- 
tween a STM superconducting tip and a bulk supercon- 
ductor, or a quantum point contact in break junctions 
[15| . The same model can also describe a quite different 
system: a suspended carbon nanotube between super- 
conductors in the Fabry-Perot regime for a transparent 
barrier [2Q|. Different groups have realized suspended 
carbon nanotubes with good mechanical properties |8|, |9| . 
Also transparent contacts between superconductors and 
nanotubes forming SQUIDs have been observed |13, [231 . 
The motion of the nanotube induces a modulation of the 
gate potential seen by the nanotube and thus modulates 
the transparency of the electronic mode. Note also that 
the mechanical coupling considered here is different from 
the one investigated in Ref. [28[ for a suspended carbon 
nanotube in presence of magnetic field. 

We will assume that the phase difference between the 
left (L) and right (R) superconducting leads, 4) = <i>L — 
(f)ji is time independent, i.e. there is no voltage drop at 
the junction. We consider the case of a single-channel 
superconducting point contact which can be described, 
for instance, by the following tight-binding Hamiltonian 
j29| (we set the units H = e = ks = I) 



H{t)= J2 Hx + VEAlit) 



(2) 



X=L.R 



Here Hx=l,r are the Hamiltonians of the isolated X 



L, R leads given by 

i£X 



E 

<i,j>ex 



{4vo^j + H.c} , (3) 



where vq is the hopping matrix between next-nearest- 
neighbor sites within the isolated X = L,R supercon- 
ductors. We use the standard notation for the Pauli ma- 
trices in Nambu space {d'x,d'y,az}- The second term in 
the r.h.s. of Eq. 1^ describes the (time dependent) hop- 
ping between the leads defined as 



VEAlit) 

v{t) 






- i^]ivHt)iJ;i 



J4>a^/2 



(4) 
(5) 



Notice that by a standard choice of the gauge, we have 
included the superconducting phase difference into the 
hopping term v. 

Wc focus our study on the electromechanical (EM) 
properties of the junction. We assume that one of the 
leads is vibrating at a radial frequency uiq, thus modulat- 
ing the hopping term between the left and the right lead. 
In the limit where the amplitude of the oscillations is 
small, the time-dependent hopping will be linearly mod- 
ulated as 



v{t) = w{l + acos(ajot)} , 



(6) 



where a = a{l/v)dv/dx <^ 1, with x the displacement 
of the oscillating lead and a the amplitude of oscillation. 
We will give estimations for a in Sec. IIVI 

For large amplitude oscillation a detailed microscopic 
model of electron transport is needed. The standard tun- 
neling picture gives a simple exponential dependence of 
the hopping on the displacement v{x) — v{xo)e^^'^^^"''''^ 
where A is the tunneling length. But this picture would 
be different for the suspended nanotube, where a linear 
dependence on the displacement is supposed to hold for 
large oscillation amplitudes. For this reason in this paper 
we concentrate on the linear displacement model, which 
demands only a single parameter ( a ) to describe the 
EM coupling. 

In the present model, the vibrational state of the junc- 
tion is considered as an external time-dependent pertur- 
bation, without its own dynamics [36| . Before proceeding 
to determine the current from the microscopic Hamilto- 
nian [Ec^. ([2)], wc start the next section by deriving an 
expression for the Josephson current within the frame- 
work of an effective two-level Hamiltonian. Wc will see 
that in a certain range of parameters, such a model pro- 
vides an accurate description of the electronic dynamics 
of the contact. 



III. DC JOSEPHSON CURRENT 

A. Andreev Two-Level Model 

In a single-channel superconducting junction, as the 
one described by Eq. ([2|), the equilibrium spectral density 



is characterized by the two ABSs |-|-) and |— ) with ener- 
gies given in Eq. ([T]). In terms of the hopping [Eq. ([5])], 
the transmission factor r of the junction is defined as 



4/3 



(1 



(7) 



where /3 = (v/vq)'^ is the ratio between the tunnel hop- 
ping amplitude and the electrode bandwidth vq. In the 
equilibrium case, the dc current is carried exclusively by 
the ABSs and can be written as the sum of two opposite 
contributions 211 



I'^l - I-n- 



I+n 



+ "-+ : 



I- = -1+ = -2—-UJA = 



A2rsin((/)) 
2ujA 



(8) 
(9) 



where n± is the occupation of the |±) ABSs which is 
given by the Fermi distribution function /(iw^i). Thus, 
in the finite temperature T case one finally obtains the 
well known expression for the Josephson current (see Ap- 
pendix) 



I'^lW 



A' 



^^^^^tanh(^). 
2uiA 2T' 



(10) 



In the limit of zero temperature, only the negative ABS 
(|— )) is populated, and contributes positively to the cur- 
rent in Eq. (|10|) . At finite temperature T, the positive 
ABS (1+)) gets populated due to the thermal smearing 
of the Fermi distribution and, according to Eq. (|8]) con- 
tributes negatively to the Josephson current. 

Now let us consider the perturbation originated by 
the mechanical oscillations. If the frequency and am- 
plitude of the perturbation are sufficiently small, one can 
still describe the current as the contribution of the two 
ABSs. In this case it is convenient to work with an ef- 
fective two-level model Hamiltonian similar to the one 
derived in Refs. [22, |30| from the microscopic Hamilto- 
nian Eq. ^. The main difference between the prob- 
lem at hand and that considered in Refs. [22, [SO] is that 
the time-dependent parameter is not the superconduct- 
ing phase, but the transparency of the junction. Adapt- 
ing the method to our problem, it gives the following 
effective time-dependent Hamiltonian for dr/dt <C A/fi: 



h{t) = A cos — CTz -f- A^/f — T{t) sin —dy 



(11) 



This Hamiltonian is written in the ballistic basis of right 
and left moving electrons that is the eigcnbasis in the 
perfectly transmitting case (t = 1). However, it is 
more convenient to write the two-level Hamiltonian in 
the instantaneous Andreev basis (23]. For that sake, 
one performs a time-dependent unitary transformation 
H = WhU -iWdU/dt, where U{t) = e^'^'Te^'^^^^^y 
and e(t) = (1/2) arctan[^l - T{t) tan 0/2] . The two- 
level Hamiltonian in the instantaneous Andreev basis is 



then given by 



Ha 



IdT 



UJAf^z 



8 dt y/l-T{t) 1 - T{t) sin2(0/2) 



(12) 

The off-diagonal terms describe the couphng between the 
Andreev levels due to the EM oscillations. The current 
operator written in the same basis reads 



Ia 






90 



ijJA 



(13) 



We consider here the time-dependent transmission deter- 
mined by expressions Eqs. (jH]) and ([7]) which in a linear 
approximation with respect to the amplitude a is 



^(t) w T + 2t\/\ — TCI cos(a;ot) . 



(14) 



In order to avoid unphysical values of r larger than one 
one has to impose a < \/l — t/2t. More over in the par- 
ticularly interesting case of very transparent channel, the 
linear term of Eq. (J14p vanishes and it would be necessary 
to consider the quadratic one. Thus, for (r — > 1) com- 
parison between the linear and quadratic term imposes 
the condition a ^ 2-^/1 — t. 

We can now apply the method developed in Ref. [23 
to obtain time-averaged quantities like the current and 
the level population close to the first resonance, i.e. 
ujq w 2uia- Within the rotating- wave approximation and 
for one-phonon assisted processes we obtain for the dc 
current 
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DC 
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(0) 
DC 
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02 



{2lja - wo)2 + ^l 



(15) 



where the Rabi frequency 

Oi^ = aTsin(0)a;oA2/(2wA)^ 



is proportional to the coupling strength a. Notice that 
at the resonance, the dc current vanishes. This is due 
to the fact that when ujq = 2uja, resonant transitions 
between the Andreev levels take place and both levels will 
be on average equally populated. According to Eq. ([5]), 
this leads to a decrease of the dc current, which vanishes 
exactly at the resonance. The width of the resonance as 
given by Eq. P^ is proportional to the amplitude a of 
the oscillation. Thus, in principle by measuring the CPR 
one could determine both the amplitude and frequency 
of the contact oscillation as given by Eqs. ([T5|) and ([T6|. 
Finally, the time averaged population of the upper and 
lower Andreev states are given by 



1 



1 



02 



n> = - 



2 {2ujA - cjo)^ + ^R 



2{2lja 



wo)2 + n\ 



(17) 



(18) 



These results are particularly simple and transparent, 
but they rely on different approximations. In the fol- 
lowing section we will thus perform a fully microscopic 



calculation in order to check their validity. The method 
used does not need a hypothesis on the slowness of the 
mechanical frequency nor on the amplitude of the oscilla- 
tions (a) and keeps the full description of the electronic 
system. The approach allows us to take into account the 
effect of the continuum spectrum and in particular will 
remove the limitation on a <C yl — r. We will see that 
the largest deviations arc exactly where higher orders of 



a become important {a > \/T 
nant (t> is close to tt. 



r) and when the reso- 



B. Nambu-Keldysh Green function method 

In the preceding section we have determined the dc 
Josephson current from an effective two-level model. In 
this section we introduce a numerical method that al- 
lows us to compute exactly the dc Josephson current in 
the presence of EM coupling from the microscopic model 
introduced in Sec. |lTl From our results, we will able to 
verify the range of validity of Eqs. ([T5|) and ([T5)) . and 
to obtain the current out of that range. The method is 
based on the computation of the Nambu-Keldysh Green 
functions (GFs) for the fermionic fields of the Hamilto- 
nian ([2]). These are defined as 



fx'it^t') = -^{Tci^xmlc'it')), (19) 



G 



where Tc means time ordering along the Keldysh contour, 
a, /? = ± denote the Keldysh branches, and X, X' = L,R 
stand for electrode indexes. From charge conservation at 
the L-R interface, one obtains the expression of the mean 
current crossing the junction at time t in terms of the GFs 



(16) /(t) = tr{a4i}(t)G+-(f,<)-f)t(f)G+-(t,i)]}. (20) 



It is convenient to rewrite the Hamiltonian [Eq. ([2])] as 
the sum of an unperturbed part plus the time-dependent 
perturbation 



where 



H = Ho + V{t) , 


(21) 


Ho = y^ Hx+vWt, 


(22) 


X = L.R 




/{t) ~ av cos{ujot)WT , 


(23) 



and 



Wt = ijU^e'^^'/^i'R + i^U.e-^^^^/^V'z 



(24) 



The corresponding Dyson equation is shown diagram- 
matically in Fig. [51 The thin lines represent the unper- 
turbed Green's functions, i.e., those associated to Hq, 
while the thick lines are the exact GFs which take into 
account the external time-dependent perturbation repre- 
sented by a wavy line. Due to the time periodicity of 



IV. DISCUSSION OF THE RESULTS 
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FIG. 2: Schematic representation of the Dyson equation 
[Eq. (|29p ] for the Green function in the Electrode-Nambu- 
Keldysh-Floquet space. The thick (thin) hues represent the 
Green functions in the vibrating (nonvibrating) case, while 
the wavy hne represents the external time-dependent pertur- 
bation. 



Following the procedure described in the previous sec- 
tion we present in Fig. [3] the numerical results obtained 
for the dc Josephson current in the presence of the EM 
interaction for a highly transmitting junction (r — 0.98). 

The dark dotted line represents the dc current Ijjq{4>) 
obtained in absence of vibrations {a ~ 0) as it is given 
by Eq. ^. For cjq = 0.4A and a = 0.04 the dashed 
and solid lines show the results of the two-level model 
[Eq. (jTSJ) ] and the full numerical solution, respectively. 
As anticipated before, the EM coupling induces an an- 



thc Hamiltonian one can write the GFs in Floquet rep- 
resentation and the corresponding current operator as a 
Fourier series 



;e— (-*')Gf^,„(..) 



/(i) = J2 e~™'^°*/n • 



(25) 
(26) 



We are interested in the dc component (n = 0) of the 
mean current given by 

Idc{<P) = 

(27) 

where the hopping matrix v in the Floquet space is de- 
fined as 



WmO 



v{h: 



m,0 



^<5m,Tl} 



(28) 



The Dyson equation in the frequency representation is 
then given by (cf. Fig. [2]) 
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Here the self-energy term is given by 



XiXv.rn 



0'A^)'"x{X^°ra■,±l 



(29) 



(30) 



and the expressions for the free propagators g'^-^,{ijj) in 
the absence of vibrations (a = 0) are given in the Ap- 
pendix. 

In the case a 7^ 0, however, the GFs G'^^,.^{(jj) can 
only be found numerically by solving the Dyson equa- 
tion (|29|) . The latter constitutes a linear system in 
the Electrode-Nambu-Keldysh-Floquet space of dimen- 
sion 2 X 2 X 2 X {2Nph + 1) and is solved by exact nu- 
merical inversion. The maximum number of vibrational 
quanta in the system N„}, is increased until convergence 
of the solution is found [33] . 




FIG. 3: Josephson current for ujq = 0.4A, r = 0.98, a = 0.04, 
T — K, and 77 = O.OOIA (77 is the inverse of the relaxation 
time of the superconductor as defined in the Appendix) . Plain 
(dashed) curves correspond to the full numerical (analytical) 
calculations. Dotted curve: current in the case of vanishing 
a. Inset: Same plot on a larger scale. 

tiresonance on the dc Josephson current when the con- 
dition luq ~ 2u)A is fulfilled, namely when the phase dif- 
ference between L and R superconductors reaches the 
critical value 



1 



Wo 



2arcsin,/-[l-(^)2] 



(31) 



The signature of such a modulation is the presence of dips 
in the current-phase relation, the position of which pro- 
vides a measure of the vibrational frequency wq through 
the resonance condition Eq. (PT|) . The width of the 
dip is proportional to the coupling a between the nano- 
resonator and the vibrational mode of the junction be- 
ing excited. It is well approximated by the expression 
(|16p. The numerics show some additional dips in the 
CPR that are associated with higher-order transitions. 
Such processes which are obviously absent from Eq. ([15]) 
could be incorporated by extending the analytical cal- 
culation to the next leading orders in powers of a (see 
Refs.|24l. [23). According to the upper panel of Fig. HI 
the resonance dip is accurately approximated both in 



position and width by the Eq. ((T5|) . The discrepancy 
between the simpUfied two-level model and the full nu- 
merical calculation becomes visible when increasing the 
value of a. This deviation is due to the fact that for 
large enough values of a the condition a < \/\ — t is 
not satisfied. In the particular case of the upper panel 
of Fig. H a = 0.01 - 0.025, r = 0.999 and wq = O.lA. 
Therefore the analytical result is valid as far as a < 0.03. 
For larger values of a, Eq. (jlSp is no longer valid and one 
has to resort to numerical results. These are shown in 
the lower panel of Fig. 21 where a = 0.04 — 0.1. The pa- 
rameters used in Fig. 4 correspond to a high-frequency 
oscillator wq/Stt = 1 GHz and a superconductor with a 
small superconducting gap A « 0.04 mcV. Suspended 
carbon nanotubes between two superconductors seem to 
be good candidates to reach this regime. For instance, 
the resonant frequency for the fundamental mode of such 
a vibrating nanotube was reported in Ref. [3l| to be in 
the range of 500 MHz. In Ref. [33 a carbon nanotube 
was connected to a superconducting Al/Pd bilayer elec- 
trode with a BCS gap of 0.08 meV which is in the range 
of our estimation. 
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FIG. 4: Upper panel: Josephson current for the cases of 
EM coupling a = 0.01, 0.0175, 0.025 (within the range of va- 
lidity of the analytic model). Plain (dashed) curves corre- 
spond to the full numerical (analytical) calculation. The ar- 
row is denoting the direction of increasing values of a. Lower 
panel: Josephson current for the cases of larger EM coupling 
a — 0.04,0.06,0.08,0.1, as obtained from the full numerical 
calculation. Common to all curves: t = 0.999, ojq = O.lA, 
rj = O.OOIA and T = K. Topmost curve: current in the case 
a = 0.0. 



As mentioned above some differences between the nu- 
merical results and the analytical one emerge for close 
to TT that become more pronounced when increasing the 
value of the coupling strength a, as shown in the lower 
panel of Fig. HI for which a = 0.04-0.1. Notice that 
in this regime the exact numerical calculation predicts a 
change of the current sign close to (f> ~ n due to higher 
order processes. 



The results presented show that measuring the current 
phase relation in a Josephson junction coupled to a me- 
chanical oscillator can allow the detection of its periodic 
oscillation (for example when the oscillator is driven by 
an external force) . We found that the CPR displays sharp 
dips when the resonance condition ujq = 2uja is met. Ac- 
cording to Eq. ([1]) for t ^- 1 and LdA{4' = tt) — >■ 0, one 
can in principle always satisfy the resonant condition. 
In reality this can be difficult, since it requires r being 
very near to one and A for standard superconductors is 
much larger than the typical mechanical frequency. This 
problem can be solved by using superconducting alloys 
with smaller gap, or by introducing a magnetic field in 
order to reduce A to a value which is slightly larger than 
the mechanical resonance. Then the fine tuning of the 
resonance with the phase bias can be possible. 

The second crucial parameter in order to observe this 
effect is the coupling constant a. A reasonable esti- 
mate of the order of magnitude of a can be obtained 
by considering the experiment of Ref. [lO|, which was 
performed on a driven nanomechanical oscillator in its 
normal metallic state. The authors of Ref. [l0| esti- 
mate the resistance dependence on the displacement to 
be {1/ R)dR/dx ^ O.lnm"^, with a typical displacement 
in the driven case of the order of a nm. Thus a very crude 
estimate of the order of magnitude of a is 0.1, that is a 
quite strong coupling. Remarkably, for a suspended car- 
bon nanotube in the Fabry-Perot regime we find a similar 
order of magnitude. This can be estimated by using the 
responsivity of the transparency to an external change 
of gate voltage from Ref. [2y| and {l/Cg)dCg/dx « l/d, 
where Cg is the gate capacitance, and d is the distance of 
the nanotube from the gate. It is clear that at this level 
these are only crude estimates of the order of magnitude 
of a, but the result is encouraging. 

In principle one could also detect thermal motion of 
the mechanical oscillator with this method if the qual- 
ity factor (Q) of the mechanical oscillator is sufficiently 
large As a matter of facts, the thermal motion can be 
seen as a sequence of periodic oscillations with a coher- 
ence time given by Q/wq. Averaging the current over a 
time much longer than this time corresponds to averag- 
ing the current obtained above over different values of 
the amplitudes of oscillation (and thus of a). One thus 
expects also in this case a dip in the CPR, but with a 
width that is controlled by the average of a. In the case 
of Ref. [10[ this motion is tiny and gives that on average 
a ~ 10~^. This gives an extremely sharp dip, and thus 
its observation is subject to a very accurate detection of 
the CPR. 

All the results presented so far are for the zero tem- 
perature limit. In the case of finite temperature, the 
coupling of the quantum point contact to the mechanical 
oscillator may lead to enhancement of the supercurrent as 
was discussed in the context of a microwave field [2J, [23 ■ 
The enhancement of the current is due to inelastic pro- 
cesses which promote particles from the continuum spec- 
trum (a; < —A) to the lower ABS. This phenomenon (not 
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FIG. 5: Josephson current in the finite temperature case 
T = 0.0, 0.5, 1.0, 2.0a;o, as obtained from the full numerical 
calculation. Parameters of the plot: r = 0.999, ljo ~ O.IA, 
a = 0.01 and tj = O.OOIA. Topmost curve: current in the 
case a — 0.0. 



mechanical oscillation [3J, |35| since the interesting region 
for this device is the very transparent case. 



V. CONCLUSIONS 

In conclusion wc have shown the possibility of detect- 
ing ultrafast oscillations of a nano-scale Josephson junc- 
tion by analyzing its dc current-phase characteristics. In 
the high transmission regime r « 1, the onset of elec- 
tromechanical coupling results in the appearance of dips 
in the Idc{<I>) characteristics for precise values of the 
phase difference 4>. The location of those dips provides 
a new way to measure the vibrational frequency wq of 
the oscillator, and their width is directly proportional to 
the EM coupling strength a. If the latter is sufficiently 
small, we have derived an effective two-level Hamiltonian 
[Eq. ([T^ ] which describes quite accurately the dynamics 
of the contact in the presence of an EM perturbation. 
Our results provide a new way of characterizing the mo- 
tion at the nanoscale. 



shown here) is an analog to the superconducting stimu- 
lation by acoustic waves in bulk materials discussed by 
Ehashberg and Ivlev in 1986 [H. 

But more important to the efficiency of this device as 
a detector is the fact that thermal fluctuations can pre- 
vent the observation of the mechanical oscillations. For 
temperatures of the order of and larger than ujq = 2uja, 
the thermal occupation of the two-level system tends to 
the value 1/2: n± = f{±ujo/2) « l/2TuJo/iST) + o{uj^). 
This effect will reduce the overall value of the Josephson 
current; nevertheless by the exact numerical solution, we 
find that a signal is still visible till moderately high tem- 
peratures. As it is shown in Fig. [5j for T = 0.5 — 2wo the 
CPRs maintain a local minimum at the resonance. The 
position of the dip is independent of the temperature [cf. 
Eq. pi[) ] while it widens for higher temperatures. 

One should emphasize that the back action of the 
Josephson junction on the mechanical oscillator is ne- 
glected in the present work. The establishment of an 
effective two-level model for the Josephson junction de- 
scription opens the way to considering the full dynam- 
ics of the two systems coupled. At this stage, we can 
evaluate the ratio of the average of the back-action 
force {Ft,a{t)) = —{dh{t)/dx) to the elastic force F^i = 
-mwgX. We find that {Fba{t))/Fei sa 3A/4ma;gA^, where 
we introduced the characteristic length A = v/{dv/dx). 
Interestingly, this ratio is independent of the amplitude 
of the oscillations. For a single wall carbon nanotube 
of mass m « 10~^^ kg that oscillates at the frequency 
wo/27r ~ 1 GHz, we rou ghly estimate this ratio to be in 



the range 10 — 10 ^ [38j. Although very small, this 



back action may lead to interesting effects as the cooling 
effect of the mechanical degree of freedom in the presence 
of an external magnetic field [231 • Note also that we have 
safely neglected the effects of Coulomb blockade on the 
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Appendix : Free Green functions 

Here we determine the GFs of the nonvibrating Joseph- 
son junction (a = 0). We first consider the isolated Right 
electrode described by the time- independent Hamiltonian 
(Eq. ((31) ). The retarded (R) and advanced (A) surface 
Green functions of such a system can be found analyti- 
cally by making use of the periodicity of the Hamiltonian 
when writing its Dyson equation 



9T''^\u.) 



Aax 



r,=i?.(A) 

^0ff_R 



{uj)4} ' (32) 



where w,; = cj + iij. The small imaginary part describes 
inelastic scattering in the leads within the relaxation time 
approximation and is the smallest energy scale of the 
problem. In the small gap limit (A,aj << 2|wo|, va is the 
electrode bandwidth), one finds for the solution of Eg 15^ 
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--RiA) 



(uj) = N,,{uj)luj,-i + Aa^i 



(33) 



where N,{c.) = ^±[^m^^ + ,r,sign{cv) ^^^^} . e{\A\-\u;\)^\A\^ - 0J|~^r^sign{u)e{\uHA\)^u?, ~ \A\ 



We connect now the R-lead to the L-lead through the 
time- independent part of the tunnel Haniiltonian (Eq|3|) . 
The Dyson equations for the retarded (advanced) GFs of 
the entire nano-junction read 



AP 
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{34) 



(35) 



Eas. p4j|55)) can be solved analytically and provide the 
expressions for the unperturbed GFs used in Eq. (|29p 
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(36) 



a{4>) -6(-0) 
6(0) ~a{~(b) 



a(0) = e'^''uj;~e-"^''ui^A (37) 

6(0) = uj^{e"t'^^A-e-'^/'^ui^) (38) 

where we introduced the intermediate function g,-i{i-u) = 



and the complex number w^ 



A 



l+/3e" 



the mod- 



lA - ^ 1+/3 

ule of which is the Andreev bound state energy 

ujA ~ AWl — rsin (0/2). The remaining components 
of the GFs are obtained by changing simultaneously 
left and right electrode indexes and the sign of the 
superconducting phase difference, i.e., (i, 0) — )> (-R, — 0). 

Finally, the non diagonal component of the Keldysh GFs 
can be found using the relation (valid at equilibrium only) 



G+-^,{u;) = f{u;){Gi^,ico) Gf^,(c.)} 



(39) 



where / is the Fermi distribution function. From this 
equation one can compute the equilibrium Josephson cur- 
rent (Eq llOp by integrating the following expression 
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n I 



dw 
2 
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} ■ (40) 
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